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We consider sterile neutrinos with rest masses ~ 0.2 GeV and with vacuum flavor mixing angles 
sin 8 > 10 for mixing with r-neutrinos, or 10 < 6 < 10 for mixing with muon neutrinos. 
Such sterile neutrinos could augment core collapse supernova shock energies by enhancing energy 
transport from the core to the vicinity of the shock front. The decay of these neutrinos could 
produce a flux of very energetic active neutrinos, detectable by future neutrino observations from 
galactic supernova. The relevant range of sterile neutrino masses and mixing angles can be probed 
in future laboratory experiments. 
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I. INTRODUCTION 

Neutrino masses are usually incorporated into 
the Standard Model (SM) by the addition of 
SU(3)xSU(2)xU(l) singlet fermions, often called right- 
handed neutrinos [l| . These gauge singlets can have Ma- 
jor ana masses as low as a few eV 0, or as large as the 
Grand Unified scale If the Majorana mass terms are 
large, the particles associated with the singlet fields are 
very heavy. However, if the Majorana masses are below 
the electroweak scale, the corresponding degrees of free- 
dom appear in the low-energy effective theory as so called 
sterile neutrinos. Light sterile neutrinos with masses of 
a few keV could be the cosmological dark matter 0, 0| , 
their production in a supernova could result in large pul- 
sar kicks IfJ , and they could affect supernovae in a variety 
of ways [6(. The same particles can play an important 
role in the formation of the first stars H , baryogenesis [1| , 
and other astrophysical phenomena [9(. 

In this paper we investigate the effect of heavier sterile 
neutrinos in astrophysical systems. We show that ster- 
ile neutrinos with masses ~ 0.2 GeV and small mixing 
sin 2 9 ~ 10~ 8 , with either the muon or tau neutrinos, 
could be produced in supernova cores and subsequently 
augment, via their decay, the energy transport from the 
core to the region around the stalled shock, thereby in- 
creasing the prospects for a core collapse supernova ex- 
plosion. This scenario could be testable, both by future 
laboratory searches for heavy sterile neutrinos and by 
observations of the neutrino signal from supernovae. 

Our analysis differs from earlier work. Heavy neutral 
lcptons, with masses 10 keV - 10 MeV, produced in su- 
pernovae, have been considered for setting limits and 
for powering supernova explosions [10|. However, the 
mass range we consider here, 145 - 250 MeV, is quali- 
tatively different in that the sterile neutrinos decay pre- 
dominantly into a pion and a light fermion. If the mixing 
angle with the electron neutrino is negligible, and the 
sterile neutrino mass M s is in the range m n o < M s < 
(m^o +to m ), the daughter pion is the neutral pion, which 
decays into two photons: — > i/^ir — > [ll[, 
where a = (fi, r). This decay mode, with lifetime ~ 0.1 s, 



changes the impact of sterile neutrinos on the supernova 
explosion. To distinguish sterile neutrinos that decay 
mainly into photons from the other types, we call them 
eosphoric. 1 As a consequence of small mixing angles, 
sin 2 9 ~ 10 ~ 8 — 10~ 7 , the eosphoric sterile neutrino pro- 
duction and propagation history inside the supernova en- 
vironment could be significantly different from those of 
the neutral leptons with electroweak scale interactions 
which have been considered previously [To| . 

We assume that the particle physics lagrangian at low 
energies comprises the Standard Model, albeit modified 
as indicated. The Standard Model was originally formu- 
lated with massless neutrinos i>i transforming as compo- 
nents of the electroweak SU(2) doublets L a (a = 1, 2, 3), 
but here we will extend it to include seesaw mass terms 
for neutrinos [l|, in which we allow the Majorana masses 
to be below the electroweak scale: 

C = C SM +iN a fiN a -y aa Hi L a N a -^- N c a N a +h.c. (1) 

The neutrino mass eigenstates 

{fx V2 v 3 > v i > v 2 1 v n } a re linear combina- 
tions of the weak eigenstates {v a , N a }. The states vf \ 3 

(sj 

are active and have masses below 0.2 eV, while v\ n 
are sterile. In particular, several recent studies focus 
on the ^MSM [4J, a model with n = 3 sterile neutrinos: 
one with a few keV rest mass (dark matter), and two 
nearly degenerate, heavier states. This model facilitates 
leptogenesis via neutrino oscillations [|[. 

We will consider sterile neutrinos with masses 145- 
250 MeV, and with vacuum flavor mixing angle sin 2 9 ~ 
10~ 8 — 10 -7 for mixing with either muon or tau neutrinos 
in vacuum. The current bounds [12J allow this range of 
mixing angles with both muon and tau neutrinos. For 
sterile neutrinos mixing only with tau neutrinos an even 
broader range is allowed [r2l . [l3| . Sterile neutrinos with 
somewhat smaller vacuum mixing could produce a pre- 
nucleosynthesis matter-dominated epoch in the early uni- 
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verse, and this is being investigated separately [14[. In 
section [TTJ we investigate the production of these neutri- 
nos in the core of a hot proto-neutron star. In section HTT1 
we discuss how the decay of eosphoric neutrinos leads to 
energy deposition in the mantle above the post-collapse 
neutron star which could help in shock revival. Section 
IIVI is devoted to the neutrino signal produced by these 
decays, and the associated observational signature of this 
scenario. We show that this is consistent with the neu- 
trino detection from SN1987A. 



II. STERILE NEUTRINO PRODUCTION 

Sterile neutrinos are produced in a hot proto-neutron 
star by electron-positron or neutrino pair annihilation, 
and by inelastic scattering of ^( T ) an d ^Vi(r) on an y °f 
the fermionic species. Since n,p, e~ , e + , v ei v e are degen- 
erate in the core, fermi-blocking will render pair annihi- 
lation of and inelastic scattering on the non-degenerate 
neutrino species v^, v T the dominant process for N s pro- 
duction [151 ]. The amplitude for these processes, in terms 
of the Mandelstam variables, is [lq : 



Process 


Amplitude 




+ - 


-» Vp(T) 


+ Vs 


sin 2 6 


32G 2 F u{u - 


M?) 


5 rW 


+ »rM ~ 


"> ^(r) 


+ Vs 


sin 2 6 


> 8G% u(u - 


M 2 ) 




+ V »{r) - 


-» ^r( M ) 


+ Vs 


sin 2 


9 8G% s{s - 


Mt) 






-» ^r( M ) 


+ v s 


sin 2 I 


> 8G% u{u - 


M 2 ) 



TABLE I: Dominant processes contributing to sterile neutrino 
production in the supernova core, for mixing with the muon 
(tau) neutrino by angle 6. Since the sterile neutrino is its own 
antiparticle, the charge-conjugated counterparts of the above 
processes also contribute to u a production. 

We have calculated the sterile neutrino production rate 
numerically by appropriate phase-space integrations in- 
volving particle distributions and blocking factors [T^ |. 
Sterile neutrinos with masses M s w 145 — 250 MeV and 
vacuum mixing sin 2 6 ~ 1CP 8 — 10~ 7 are not trapped; 
they stream freely out of the supernova core. Since both 
the production cross sections a ~ sin 2 G 2 F E 2 tx T 2 and 
the ffj.iT) number density n ~ T 3 grow rapidly with tem- 
perature, the sterile neutrino flux and luminosity depend 
strongly on temperature. Our numerical results, for a 
sterile neutrino with mass M s = 200 MeV, and for tem- 
peratures in the range 15-70 MeV, are well approximated 
by the following fitting functions: 
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FIG. 1: Energy spectrum of sterile neutrinos, of mass M 3 = 
200 MeV and mixing angle sin 2 6 — 5 x 10~ 8 ; produced in the 
temperature profiles shown in Fig. 9 of Ref. [lq ]. 
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TABLE II: Total luminosity in sterile neutrinos, for the tem- 
perature profiles of model GM3 of Ref. The mixing angle 
is taken to be sin 2 8 = 5 x 10 _s . 



Supernova cores typically reach peak temperatures 
T ~ 35 MeV, at R ~ 10 - 20 km, for a time interval 
of a few seconds [13] ■ Equations ^ and ([3]) then imply 
that the sterile neutrino burst will last 1-5 s and drain 
total energy E s ~ 10 51 erg. Sterile neutrinos will be 
emitted with average Lorentz factor 
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which is rather insensitive to temperature within the rel- 
evant range. 

In particular, we employ models of supernova cores 
from 1-dimensional simulations: the model GM3np of 
a proto-neutron star, from Ref. and the model 

sl5Gio_ld.b, from Ref. [lj| , for a star with main-sequence 
mass 15 Mq. This simulation implements general- 
relativistic gravity advanced description of neutrino in- 
teractions and full spectral treatment of neutrino trans- 
port, but it fails to generate a successful explosion. We 
reproduce the temperature profiles of these models and 
present the resulting sterile neutrino luminosity-energy 
spectra in figures [1] and [21 Tables [Til Mil contain our nu- 
merical results for the total luminosity in heavy sterile 
neutrinos. 

These calculations show that sterile neutrinos in the 
parameter range under consideration can play an impor- 
tant role in the energy transport and deposition budget 
in the supernova environment. Although eosphoric ster- 
ile neutrino production is suppressed by their small mix- 
ing angles and relatively large rest masses, their poten- 
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FIG. 2: Energy spectrum of sterile neutrinos, of mass M 3 — 
200 MeV and mixing angle sin 2 8 = 5 x 10 -8 , produced in 
model sl5Gio_ld.b of Ref. -IS] (shown in Fig.20 of Ref. QJ]) 
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TABLE III: Total luminosity in sterile neutrinos, for the tem- 
perature profiles of sl5Gio_ld.b of Ref. [lj|. The mixing angle 
is taken to be sin 2 8 = 5 x 10 -8 . 



tial impact is enhanced by their free-streaming out from 
the core. The net result, for typical supernova cores, 
is a total energy output in sterile neutrinos of the or- 
der of the initial bounce-shock energy, ~ 10 51 erg. As 
a result, eosphoric sterile neutrinos carry an amount of 
energy that, if deposited under or in the vicinity of the 
shock, could have an impact on the supernova explosion 
mechanism. 



III. SUPERNOVA SHOCK ENHANCEMENT 

In standard gravitational collapse, a shock forms at 
the edge of a cold (T » 2 MeV), low entropy, homol- 
ogous "inner" core piston and then subsequently propa- 
gates outward through the outer core, heating this region, 
but losing energy through photo-dissociation of heavy 
nuclei into free nucleons and alpha particles and, as a 
consequence, in ~ 100 ms evolves into a "stalled" stand- 
ing accretion shock. With standard weak interaction 
physics, this shock may be revived by neutrino heating 
or neutrino heating aided by multi-dimensional hydrody- 
namic effects, ultimately producing a supernova explo- 
sion [llj [2(| . The inner core heats up and the whole core 
de-leptonizes on a time scale ~ 1 s, the neutrino diffusion 
time scale. 

The eosphoric neutrinos decay into 7-ray photons and 
very energetic (~ 80 MeV) active neutrinos outside the 
neutron star but inside the envelope. These decays in- 
crease the total energy of the envelope, making it easier 



to eject the material from the gravitational potential well 
of the neutron star. Their decay time is [llT ] 
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where f v = 131 MeV. For mixing with v^ 7 in the window 
M s = 145 - 250 MeV and sin 2 9 = 10~ 8 - 10~ 7 , their 
decay time is r s ^0.1 s. Emitted from the core with 7 S ~ 
1.5, they would deposit their energy at distances R s ~ 
10 9 cm. For mixing with v T , a much wider parameter 
range is allowed. A requirement that the energy output in 
sterile neutrinos remains E s ~ 10 51 — 10 52 erg, essentially 
to stay within observational supernova bounds, would 
limit the parameter region of interest, but it would still 
allow for lifetimes as small as r s ~ 10~ 3 s and R s ~ 
10 7 cm. From eqs. ([3]) and we see that this can be the 
case, for example, if M s ~ 300 MeV and sin 2 9 ~ 10~ 6 . 

The 7-ray photons quickly thermalize with matter and 
could generate a pressure gradient of sufficient magni- 
tude to affect the standard picture of the gravitational 
collapse. The energy density deposited in photons is 
m 7 = (L s /4-7ri? 2 /3 s )(l — e -R/ R =} an d the resulting pres- 
sure gradient is 
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in a region with M ~ 0.1 M Q could give an increase in 
entropy-per-baryon As ~ E s / (TM/m p ) ~ a few units of 
Boltzmann's constant per baryon. If As > 3, nuclei in 



nuclear statistical equilibrium would be melted 21[ and 
at least some of the photo-dissociation burden on the 
shock could be alleviated. Whether by helping revive 
the stalled shock or by altering the thermal environment 
in its vicinity, the prospects for an explosion likely would 
be enhanced in this scenario, even in a simplistic one- 
dimensional supernova model. 

The energy E s delivered by sterile neutrinos can vary 
greatly, depending on the core temperature, which need 
not be the same for all supernovae. Supernovae with 
E s ~ 10 51 — 10 52 erg should have a healthy shock wave 
and should be more likely to produce a neutron star, 
while supernovae with E s <C 10 51 erg are more likely to 
produce a black hole, due to a weaker shock and more 
infalling material. In any case, the existence of eospho- 
ric sterile neutrinos with the properties we consider could 
alter the standard core collapse supernova energetics con- 
siderations while staying within existing observational 
bounds. 



IV. NEUTRINO SIGNATURE 

Eosphoric sterile neutrino decay will produce or v T 
neutrinos with energies e„ ~ 80 MeV. Although very en- 
ergetic, these neutrinos likely still will escape from the 
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supernova. Their optical depth for scattering on the sur- 
rounding ambient matter is 

t v = G%el-^R < 1(T 2 , (7) 

m„ 

for p < 10 5 g/cm 3 at R ~ 10 9 cm. The optical depth 
could be somewhat greater for the case of mixing with tau 
neutrinos, because the decays of sterile neutrinos could 
take place closer to the core, as discussed above. 

Neutrino observations from supernova explosions can 
help test the scenario presented in this paper. In the 
event of a galactic supernova, experiments could de- 
tect the initial burst, 1-5 s duration, of energetic v^^ 
from eosphoric neutrino decays, followed by a longer (10- 
— 15 s) neutrino signal of ~ 15 MeV neutrinos. This 
signal will, of course, be modified by standard neutrino 
self-coupling and matter-affected oscillations 22]. The 
detection probability of the very energetic neutrinos is 
larger than that of ordinary neutrinos by a factor of 
~ (80MeV/15MeV) 2 = 28. However, they carry only 
a small portion, 10 51 erg/3 x 10 53 crg ~ 0.3%, of the to- 
tal energy carried away by neutrinos. Thus, their time- 
integrated flux will be a small fraction of the total neu- 
trino flux: 

(Jao/Ftotai) ~ 28 x 0.3% ~ 8% (8) 

The predicted neutrino signal is in agreement with the 
neutrino observations from SN1987A. The Kamiokandc 
detector observed 12 neutrino-induced events associated 
with SN1987A [23[. This implies that it could have seen 
one event originating from sterile neutrino decays. This 
is not in contradition with the data [231] . In addition, 
the emission of sterile neutrinos depends very sensitively 
on the core temperature, growing as a high power of 
this quantity. Since SN1987A was an unusual supernova, 
both from the standpoint of its progenitor star and the 
apparent absence of a pulsar in the remnant, it is also 
possible that the eosphoric neutrino production was in- 
efficient due to a lower than average temperature in the 
core. In this case an even lower flux of energetic neutri- 
nos would have been expected. We emphasize, however, 
that even with the core temperature expected in an aver- 
age supernova, SN1987A should have produced less than 
one event stemming from sterile neutrinos with the pa- 
rameter ranges considered here. The non-observation of 



gamma rays from SN1987A [24] is consistent with our 
model because, in the core collapse of a blue supergiant, 
the gamma rays should have been absorbed by the enve- 
lope. 

Modern and future detectors are expected to observe 
high neutrino-induced event rates from a galactic super- 
nova. 



V. CONCLUSIONS 

Heavy sterile neutrinos could prove to be very impor- 
tant in compact astrophysical systems, such as super- 
novae, essentially because they could move significant 
amounts of energy around in ways that ordinary active 
neutrinos and hydrodynamic motions cannot. In this 
paper, we have showed that sterile neutrinos with mass 
r~> 200 MeV and small mixing sin 2 9 ~ 10~ 8 with either 

or v T could facilitate energy transport from the super- 
nova core to the shock front, possibly ultimately leading 
to a successful explosion, but in any case altering super- 
nova energetics in ways which change the standard core 
collapse paradigm, yet produce signatures and behaviors 
that could remain within existing observational bounds. 

The eosphoric sterile neutrino scenario can be tested 
by future neutrino observations from galactic supernova 
explosions. A supernova should produce a short burst 
of very energetic neutrinos followed by a longer (and 
much more powerful) signal of lower-energy neutrinos. 
The non-observation of very energetic neutrinos from 
SN1987A is consistent with our model because of the 
low flux. However, present and future neutrino detectors 
should be well positioned to detect the (less numerous) 
80 MeV neutrinos. Observations of a Galactic supernova 
neutrino signal in these detectors could be used to con- 
strain the eosphoric sterile neutrino parameter space in 
ways which could extend or be complimentary to future 
laboratory-based neutrino mixing probes. 
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